Dephasing times in quantum dots due to elastic LO phonon-carrier collisions 
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Interpretation of experiments on quantum dot (QD) lasers presents a challenge: the phonon bottle- 
neck, which should strongly suppress relaxation and dephasing of the discrete energy states, often 
seems to be inoperative. We suggest and develop a theory for an intrinsic mechanism for dephasing 
in QD's: second-order elastic interaction between quantum dot charge carriers and LO-phonons. 
The calculated dephasing times are of the order of 200 fs at room temperature, consistent with ex- 
periments. The phonon bottleneck thus does not prevent significant room temperature dephasing. 
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Semiconductor Quantum Dots (QDs) may pave the 
way to optoelectronic devices with vastly superior per- 
formance compared to present devices [ [Q]. A detailed 
understanding of QD's optical properties is therefore of 
utmost importance, both from a basic science and from a 
technological point of view. The homogeneous linewidth 
r of an optical transition, and the dephasing time, or, 
equivalently, the polarization relaxation time, T2 = 2/T, 
are basic characteristics for interaction of light with quan- 
tum systems. In lasers, the homogeneous linewidth can 
strongly effect the laser gain. The dephasing time de- 
fines the time scale on which coherent interaction of light 
with medium takes place [ 0], and, in particular, it gives 
the ultimate time scale for realization of coherent control 
in a quantum system (see, for instance, Ref. [ B] and 
references therein). 

In bulk and Quantum Well (QW) semiconductors the 
dephasing time is usually taken to be the intraband re- 
laxation time [ Q]. At low carrier densities, this time is 
set by LO phonon-carrier scattering while at higher car- 
rier densities carrier-carrier collisions are also important. 
The energy spectrum of carriers in QDs is discrete, and 
in this sense the interaction processes in QDs are qual- 
itatively different from those in QWs or bulk materials. 
For example, if the LO phonon energy does not coin- 
cide with the separation of energy levels in a QD, LO 
phonon-carrier interaction can not lead to carrier relax- 
ation between these levels (this is the so called "phonon 
bottleneck" [ At first sight, then, in the absence 

of this efficient mechanism much longer relaxation and 
dephasing times could be expected. 

Carrier relaxation is a result of inelastic carrier-phonon 
and carrier-carrier collisions. In contrast, the dephasing 
time T2 is influenced by both inelastic and elastic col- 
lisions. Uskov et al. [ 0] demonstrated recently that 
elastic collisions of two-dimensional (2D) carriers with 
carriers confined in self-assembled QDs can lead to sub- 
stantial collisional broadening (T2 ~ 0.1 - 1 ps at moder- 
ate densities of 2D carriers) of QD spectral lines. In this 
work we study elastic collisions between QD-carriers and 



LO-phonons: these collisions can disturb the phase of 
the carrier wave functions in QDs without changing the 
populations of the carrier energy levels, and accordingly 
QD spectral lines will be broadened. We will show that 
these processes can lead to dephasing times T2 of some 
hundreds of femtoseconds for typical QDs at room tem- 
perature. These values are of the order of the dephasing 
times in bulk materials and QWs, and are in accordance 
with recent room temperature measurements [0]. 

Dephasing and broadening of spectral lines in QDs at 
low temperatures (< 50 K) are usually attributed to in- 
teractions of QD carriers with acoustic phonons [ ||s|-|l0||1. 
In this case, the continuum of acoustic phonons leads to 
a continuous band of satellites around the Zero Phonon 
Line (ZPL). This band is considered as homogeneous 
broadening of the formerly discrete QD spectral line, and 
the dephasing time is calculated from the width of the 
band. The present paper considers broadening and de- 
phasing via a polar coupling to LO phonons with fixed 
frequency. The broadening of the ZPL by LO phonons is 
interpreted as follows. QD carriers can virtually absorb 
and then emit (or emit and then absorb) LO phonons. 
This second order process does not change the final en- 
ergy of the carrier, but it does change the phase of the 
carrier wave function. This change of the phase, in turn, 
implies dephasing of the dipole for the considered optical 
transition, and accordingly a broadening of the spectral 
line [0]. 

The scattering against LO phonons gives rise to fluctu- 
ations AE(t) in the transition energy between the elec- 
tron and hole QD levels: E(t) = E<® + E s + AE(t), 
where E 1 - ^ is the unperturbed transition energy, and E s 
is the energy-shift due to the LO phonons. For simplicity, 
we ignore the Coulombic electron-hole interaction. In a 
classical treatment the line structure function of the op- 
tical transition is given by t he Fourier tra nsform of the 
correlation function ipf(r) — f(t + r)/*(i) for the oscil- 
lator f(t), (here the ovcrline indicates an average over 
fluctuations) 
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f{t)cxexp{-(i/h)[(E^ +E s )t + [ dhAEih)]} 



The correlation function ip t (r) can be expressed as 
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In this work we must use the quantum mechanical ana- 
log of Eq. (0): the fluctuations in the transition energy 
are due to interactions with optical phonons, which obey 
Bose statistics. Thus, instead of V'/( T ) we consider [ [ pi] ] 



A(t) = (f exp{-(i/h)[EW T + £ T dt'H(t')]}\ , 



where f is the time-ordering operator and H(t) is the 
effective interaction Hamiltonian describing the fluctua- 
tions in the transition energy. (• • •) indicates an average 
over the phonon ensemble. The result of our calcula- 
tions, to be described below, is that the lineshape func- 
tion, which is just the Fourier transform of (0), can be 
expressed in a compact form: 
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The shift E s and the parameter ha, which has the mean- 
ing of average quadratic fluctuation of the spectral line, 
depend on temperature, phonon interaction mechanism, 
and quantum dot geometry, and explicit expressions for 
them are given below. r LO = l/2j LO is the life-time of 
optical phonons, which is finite because of various in- 
teraction mechanisms (such as phonon-phonon scatter- 
ing or boundary scattering), and it enters our theory 
as a parameter which must be calculated separately, or 
extracted from experiment. If a <C 2 r y LO , A(uS) has a 
Lorentzian line form with FWHM T — <j 2 /^ lo , while if 
a 3> 27 i0 Eq. (|4|) gives rise to a Gaussian line with width 
r = 2\/2ln~2 cr. 

In order to arrive at Eq. (||) the following steps were 
needed. First, we derive an effective Hamiltonian describ- 
ing the interaction between quantum dot charge carriers 
and optical phonons. Next, the time-ordered expectation 
value (||) must be evaluated. Finally, the developed for- 
malism will be applied to a concrete model of a quantum 
dot, and the temperature dependence of T is computed. 

The carrier-phonon interaction is (x = e for electrons, 
x — h for holes) 



^(r,t) = ^^[a q (t) e ^ r + h.c] 



(5) 



where a q is the annihilation operator for the LO phonon 
with the wave vector q. For polar carrier-phonon cou- 
pling C| = — Cq. The finite life-time of the optical 
phonons gives rise to the following expression for the 
phonon Green function [ p3[]: 



D(q,r) = (f{[a q (r) + a q (r)][a q (0) + a q (0)]}> 



_ p — lLO\T\ 



{n + l)e 



iui L0 \r\ 



(6) 



where n = l/[ex.p(fru) LO /kBT) — 1], and u> LO is the LO 
phonon frequency. Assuming that the phonon energy 
Tilj lo is much less than the energy separation between 
the QD levels, an effective Hamiltonian H(t) can be de- 
rived by neglecting real transitions beween the levels [ S], 
and applying perturbation theory in the carrier-phonon 
interaction. Up to second order one finds 



(3) = ]T [/(q)a q (f) + /*(q)at (*)] , 
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^ (2) w = -Eh (t)aq ' ( ^ (q '" q,) 



-a q (t)a q , (t)F(q, q') + a q (i)a q , (t)-F(-q, ~q') 
-a q (t)a q ,(^(-q,q') 
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E® p (E® ip ) is the unperturbed electron (hole) energy with 
the wave function ip ep (iphp), which enters via the matrix 
element 



M xp 



dr C P ( r ) exp(iq • r)if> xl (r) 



(11) 
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where E g is the energy gap of the QD 



) is just the standard 
while the second order 



We consider broadening of the transition between elec- 
tron and hole ground states (p = 1), so that E^ = 
E g + E° el 
material. 

The first order term Eq^_(|7 
Huang-Rhys Hamiltonian 

term (||) is an extension of this theory. A similar second 
order term has earlier been shown to lead to broaden- 
of the ZPL in case of impurities in doped crystals 
iTJ]. Below we show that the Hamiltonian (^) - (||), 
which only involves virtual transitions between the QD 
levels, leads to a broadening of the ZPL in QD's in ac- 
cordance with experiments [ R]. The existence of the 
nearby levels is essential for obtaining the broadening. 
The effects of mixing of the QD levels due to transitions 
between the levels on intensities of phonon satellites have 
recently been considered in [ ITS]]. 
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The evaluation of (J3J) with the effective Hamiltonian 
(0) ~ (§)■ can be carried out exactly using the cumulant 
technique [ p5| ]. The result is (note that cross- terms 
vanish because they involve unequal number of phonons) 



A(t) =cxp{ + ^ (H^)]{T/h) (12) 

3=1,2 

f d * 2 E (T'AffW(t 1 )AF«(t 2 ))} , 
In Jo Jo ._ 12 J 

where AH^(t) = H^(t) - (H^). It is interesting to 
note the formal similarity with the classical result Eq. (||). 
The phonon averages can be carried out with the help of 
Wick's theorem, which stipulates that the various op- 
erator averages implicit in ( fl2"| ) are given as the sum 
of all pairwise averages. The time-ordered expectation 
value due to is just J^dti J* dt 2 |/(q)| 2 D(q,ii — 
t 2 ); this expression is well-known, say, from stud- 
ies of phonon-assisted tunneling in resonant-tunneling 
diodes [ Kg]]. In our case this term gives the Huang- 



Rhys shift E, 



(i) 



-Snui LO of the spectral line, where 



S = \f(q)/(nw LO )\ 2 is the Huang-Rhys parameter 
pyi3"||], as well as the satellites at energies ±nfiuj LO 
from the main line. For the present purposes it is im- 
portant to realize that the first-order term can- 
not lead to broadening of the ZPL if Jlo/^lo <C 1, 
see also Ref.[ [0]. A straightforward calculation shows 
that the contribution due to is obtained by a 

replacement |/(q)| 2 £(£i - t 2 ) -> F(q, q')[F(q', q) + 
F(— q, -q')]D(q, ti — t 2 )D(q',ti — t 2 ). This expression 
generates two kinds of terms, (i) Terms proportional 
to exp[±2iw LO (ti - t 2 )]; together with (H^) in @ 
they give the spectral shift in the second-order theory. 
They are also relevant for the two-phonon satellites at 
±n2huj LO of the main line, which however is not our fo- 
cus here, (ii) Terms without these fast time-oscillations 
which determine the zero-phonon linewidth. Collecting 
all the terms describing the spectral shift and broaden- 
ing, we find 



=-{2n + l)(S 1 + S 2 )n<j LO , 
a = \J An(n + 1)S 2 uj lo , 



where 



(13) 
(14) 

(15) 



q 



5 2 -^f(q,q')[F(q',q)+F(-q,-q')]/2(^ ) 2 . (16) 
qq' 

The second-order shift (|l^) together with the first-order 
shift E^ gives the total shift E s of the ZPL in (|). The 
parameters S\ and S 2 are extensions of the parameter 
S of the Huang-Rhys theory. Our numerical calcula- 
tions show that Si,S 2 <C 1, and also S% > S 2 . Eq. (Ill) 
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FIG. 1. The line-width as a function of temperature. The 
following parameters were used: — 10.9, e = 12.9, 

m e + ruh = 0.35, R = H — 5 nm. The solid curve is the 
result of a full numerical calculation, while the Lorentzian 
and Gaussian approximations are given as dotted and dashed 
lines, respectively. The inset shows the function g(R/H). 

together with (Q) allows an analysis of the spectral line- 
shape, and, in particular, its width as a function of tem- 
perature. 

As a practical application, we consider a cylindrical 
quantum dot of radius R and height H , with infinite con- 
fining potential. In the calculation of the average fluctu- 
ation a, we include in Eq. (10) only the leading terms in 
the sum over the states p . If 2R > H, two states have 
the same energy, and are the closest to the ground state. 
Taking into account this two-fold degeneracy, we get 



(17) 



where AE = h 2 /(2{m e + m/ l )i? 2 )[a 2 1 — a 2 ]. a lrn is the 
l-th root of the Bessel function of the order m. The 
material parameter K = e 2 (l/e OQ — l/e)/(7r 2 eo) has the 
value 26.0 meV-nm for GaAs. The dimensionless function 
g(H/R) is shown in the inset of Fig. 1. 

In the numerical calculation of the width T (see Fig. 
1) we use t lo — 1.5 ps [ [Q] so that h~f LO = 0.22 meV 
and {oj lo t LO )- J ~ 0.01. From @ 

one finds that at 

T = 300 K the average quadratic fluctuation is ha ~ 2.2 
meV (corresponding to S 2 ~ 2 ■ 10~ 3 in Eq. (fl4|)), so that 
the condition a 3> 27^0 is satisfied, and consequently 
the spectral line has a Gaussian shape. At lower temper- 
atures a decreases and can reach the regime a <§; 2-f LO , 
where the Lorentzian shape obtains. At T = 300 K the 
linewidth UT is 5.1 meV, corresponding to T 2 ~ 250 fs. 
Taking into account higher excited states in the sum (|Io| ) 
in evaluating a typically reduces T 2 by a factor of 1.5 — 2. 
Thus, one can expect room temperature dephasing times 
of the order of 100 — 200 fs. This is in good agreement 
with experimental values (200 - 300 fs) [ @]. The LO- 
phonon contribution to the linewidth is proportional to 
n(n + 1) if r <C 1/t^q, and thus vanishes as the tem- 
perature tends to zero. For example, for T < 50 K the 
calculated values of UT are below 20 /xeV (see Figure 
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1), which is below the linewidths found in experiment 
(~ 50-150 AieV) [ §0]. On the other hand, the acous- 
ticphonon related linewidth does not vanish in this limit 
[ ||], and we conclude that the LO-mechanism is dom- 
inant only at elevated temperatures, say T > 100 K. 
The present broadening mechanism also provides a pos- 
sible way of understanding fast carrier relaxation in QDs 
at elevated temperatures [ JpJ]: due to the substantial 
broading of the levels their differences no longer have to 
match strictly the LO-phonon energy, thereby weakening 
the phonon bottleneck effect. 

It is also instructive to consider the classical limit of 
our calculations. Then n» 1, and thus for T -C 1/t lo 
we find Tan 2 . This result also follows from Eq. (g) if 
we treat the phonons as classical random fields perturb- 
ing the QD carriers. This classical result allows one to 
establish an analogy with broadening in atomic gases [ 
p0|]. There, elastic atomic collisions lead to dephasing 
of optical transitions, and the homogeneous linewidth in 
the Lorentzian limit is proportional to the density of the 
colliding particles. In our case the broadening of ZPL 
is the result of elastic second order interactions between 
phonons and the carriers, and hence T (xn 2 . Note also 
that the LO phonon lifetime t lo plays the same role as 
the collisional time in atomic gases [ [p0[]. 

We finally consider the effect of the shape and size of 
the QD on the dephasing time. Eq. ( p"7| ) implies that 
a cx L, where L characterizes the size of the QD. In 
the Lorentz case V cx a 2 cx L 2 , while in the Gaussian 
case r cx a cx L. From Fig. 1 we can see that V can 
be modified by a factor 2-3 by changing the ratio H/R 
of the cylinder. Note also that a and Y depend on the 
carrier masses, a cx (m e + m^). These characteristic 
dependencies point towards a possibility to engineer the 
spectral linewidths. 

In conclusion, we have developed a theory for dephas- 
ing of an optical transition in QDs due to second-order 
elastic interaction with LO-phonons. The theory results 
in expressions which are easily evaluated numerically, and 
can also be used to draw qualitative conclusions about 
the dependence of the dephasing time on physical pa- 
rameters defining the QD. At room temperature, despite 
of an apparent phonon bottleneck, the considered mech- 
anism leads to dephasing times of the order of 200 fs, 
i.e. of the same order as in bulk, and also as observed in 
QDs. 
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